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Multivariable generalizations of the classical Hermite, Laguerre and Jacobi 
polynomials occur as the polynomial part of the eigenfunctions of certain 
Schrodinger operators for Calogero-Sutherland-type quantum systems. For 
the generalized Hermite and Laguerre polynomials the multidimensional ana- 
logues of many classical results regarding generating functions, differentiation 
and integration formulas, recurrence relations and summation theorems are 



O 
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obtained. We use this and related theory to evaluate the global limit of the 
ground state density, obtaining in the Hermite case the Wigner semi-circle law, 
and to give an explicit solution for an initial value problem in the Hermite and 
Laguerre case. 
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1 Introduction 

The Calogero-Sutherland model refers to exactly solvable quantum many body systems 
in one-dimension with pair potentials proportional to 1/r 2 (in some asymptotic limit at 
least), and which have exact BDJ-type ground states: 

N 

^o = n/i( x i) n h(xj,x k ). (i.i) 

i=i i<j<k<N 

Three particular quantum many body systems of the Calogero-Sutherland type are spec- 
ified by the Schrodinger operators 

h {h) = -^+7^f + » /2 -" £ urh^ (L2a) 

j=l U - L j j=l l<j<k<N V X J djk > 



;i.2b) 
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N • 1 i 9 
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+2/3(/3/2 - 1) E ,-T/ 3 (l-2c) 
^ (sm 6j - sm 



The superscripts (if), (L), (J) stand for Hermite, Laguerre and Jacobi respectively, and 
are chosen because of the relationship of these Schrodinger operators to generalizations 
of the corresponding classical polynomials. 

A direct calculation shows that there are eigenfunctions of the form e~^ w l 2 where 



1 N 

W (H) = ^E^?- E log\x k - Xj \ (1.3a) 

3=1 l<j<k<N 
1 TV /TV 

WW = ^E^-^Elog^ 2 - E \og\xl-x]\ (1.3b) 

a' N b' N 

W {J > = --E lo gshi 2 0j - -Elogcos 2 ^- - E log | sin 2 0j -sin 2 fc | (1.3c) 

^ i=l ^ j=l l<j<k<N 

Since these eigenfunctions are non-negative they correspond to the ground state wave- 
function ipo (i.e. they are the eigenfunctions with the most negative eigenvalue Eq). Notice 
that ipo is indeed of the type ( \L.1\) . 

Conjugation of the Schrodinger operators by the reciprocal of the ground state e^ w l 2 
gives the Fokker-Planck operators 



^-t m '^- E ^ m ^%k^-^) (L4) 

(for the coordinates Xj are to replaced by 6j). Thus the Schrodinger equation 

i^iixjht) = H^ix^t), (1.5) 
with ijj = e~* E ° t e /3W///2 P and t = r/i/3 transforms to the Fokker-Planck equation 



The Fokker-Planck equation (|1.6|) describes the evolution of a classical gas in one-dimension 
with potential energy W undergoing Brownian motion. 

Two classes of problems associated with the Schrodinger operators (|1.2|) or equivalently 
the Fokker-Planck operator (O) with W given by ( |1.3| ), are the topic of this paper. The 
first is the discussion of some mathematical properties relating to the eigenfunctions, while 
the second is the evaluation of the density in the ground state and the exact solution of 
( |1.6| ) for certain initial conditions. These problems are in fact inter-related; we find that 
the density for each system can be written in terms of a certain eigenstate and that a 
summation theorem for the eigenstates gives an exact solution of ( |1.6| ). 
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A feature of the Schrodinger operators Q1.2Q is that after conjugation with the ground 
state: 

- e ^ {H -E 0)e -^^^-^JL) (1.7) 

the resulting differential operator has a complete set of polynomial eigenfunctions. In 
Section 2 we consider the form of the expansion of these polynomials in terms of some 
different bases of symmetric functions. We note that in the N = 1 case, after a suitable 
change of variables, the operator ( |1.7| ) with W given by ( |1.3| ) is the eigenoperator for the 
classical Hermite, Laguerre and Jacobi polynomials. Previous studies of the operator for 
general iV in the Jacobi case have established an orthogonality relation. Since the 
polynomials in the Hermite and Laguerre cases are limiting cases of these generalized Ja- 
cobi polynomials, we can obtain the corresponding orthogonality relations via the limiting 
procedure. 

The generalized Hermite polynomials, which are the polynomial eigenfunctions of (|1.4j) 
with W = as given by ( |1.3a| ), are studied in Section 3. Many higher-dimensional 

analogues of properties of the classical Hermite polynomials are obtained, including a gen- 
erating function formula, differentiation and integration formulas, a summation theorem 
and recurrence relations. An analogous study of the generalized Laguerre polynomials is 
performed in Section 4. In Section 5 we relate the problem of computing the ground state 
density for the Schrodinger operators ( |1.2| ) to the computation of particular eigenstates. 
By using integral formulas for these eigenstates we are able to compute the global density 



limit for even values of the coupling (3. In the case of the Schrodinger operator ( |1.2a| ), 
the limiting global density is the well known Wigner semi-circle law. Also in Section 5, 
we give interpretation to results obtained in Sections 3 and 4 for a summation formula. 
The interpretation is in terms of the solution of an initial value problem associated with 
the Schrodinger equation ( |1.5| ) 

We conclude in Section 6 by identifying the formulas contained herein which are to 
be found in previous works, and give reference to these works (two of the most important 
references in this regard are unpublished, handwritten manuscripts). In the Appendix we 
present some results relating to generalized hypergeometric functions depending on two 
sets of variables which are of relevance to the working in Sections 3 and 4. 



2 Inter-relationships 

Let us begin by explicitly calculating the operator (|1.7| ) for W given by (|1.3|) . In all cases 



it is convenient to first change variables: for W = set yj := J p/2xj : for W = 

set yj = while for W = set yj = sin 2 <pj. We then obtain 

fj(H) ._ _^L e (3W {H) /2^ H (H) _ E )e~ pw{H)/2 
P 

",d 2 a 2 " i a > 

= VSwr^w^-^^w)- (2 ' la) 

jj{L) ;= _lpwW/2( H (L)_ Eo)e -f3 W m/2 
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E 



where 



a := {pa' - l)/2, 6 := ((3b' - l)/2, a := 2/(3. 

In the one- variable case (N = 1), these operators have a complete set of polynomial 
eigenf unctions given by the classical Hermite, Laguerre and Jacobi polynomials 

e'ft-D := (-l)»(" + a )i: ( " + °!n +1)j Q(-^ (2.2c) 

respectively, where {u) n := u(u + 1) . . . (u + n — 1). 

It is also true for general N that there is a complete set of polynomial eigenf unctions for 
each of the operators ( |2.1| ). This can be seen by computing their action on the monomial 
symmetric polynomial m K , where k denotes a partition consisting of N parts Kj. We 
obtain series of the form 

e {H) {K,a)m K + £ bfjm^ (2.3a) 
e^ L \K,a)m K + £ b^m, (2.3b) 
e iJ \K,a)m K + E + E 6 £J"V ( 2 -3c) 

respectively, where the notation < \k\ means J^jLi^j < Y^jLi K j, while the notation 
fi < k means fi ^ k but 

N N p p 

E h = E K j anc ^ E — E k j f° r eacn p = i, ■ ■ ■ , n 

3=1 j=l j=l 3=1 

Also, a^K, are coefficients independent of and 

e (fl) (s,a) = -2|/c|, e (i) (/t, a) = -\k\ (2.4) 

(the explicit value of e^(n,a) can also be computed, however it is not needed in our 
subsequent discussion). This means there are eigenfunctions of the form 

a,Wm K + ~ h fM» ( 2 - 5a ) 



4 



a&W + E 5£ } 



+ E + e &s 



(2.5b) 
(2.5c) 



M<l«l 



with eigenvalues e^(/c, at), e^ L ^(/t, a) and e( J )(/c, a) respectively. 

Rather than study the eigenfunctions in the form (|2.5|) , previous studies p4), have 
shown that it is advantageous to change basis from the monomial symmetric polynomials 
to the Jack polynomials |30], We recall that the Jack polynomial Jjf \zi, . . . , zn) is 
the unique (up to normalization) symmetric eigenfunction of the operator 



A' 



A? 



d 2 2 
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a ^ ^ - 2 fe a^- 



(2.6) 



which has an expansion of the form 



fJ,<K, 



(2.7) 



« ! m 



The notation J( a ) is usually used for the particular normalization an>h K 
However, for our purposes it is more convenient to choose a different normalization, and 
to denote the corresponding Jack polynomial by as in e.g. [ 17| . This normalization 
is specified by requiring 



( Xl +...+x N r= e ci a \ Xl ,...,x N ) 

\k\=71 

It is known (see e.g. [0) that 4 Q ) and are related by 



(2.8) 



where 



with 



Ci a \ Xl , ...,x N )= a^\ K \\j- 1 4 a \x 1 , ...,x N ) 



n k(s) k(s) 



(2.9) 



(2.10a) 



3k ■ 

K(s) := l K (s) + a(a K (s) + 1) h«(s) := l K (s) + 1 + aa K (s) (2.10b) 



In (|2.10a|) and (2.10b), k is regarded as a diagram, s denotes a node in the diagram and 
o>k(s) (Ik(s)) denotes the arm length (leg length) of the node (see e.g. [EBJ) . In terms of 
the Jack polynomials, it is known [[24|, [H], |23|, |22[| that for each partition k there is an 
eigenfunction of the form 

(2.11a) 
(2.11b) 
(2.11c) 



H K ( yi ,. 


■ ,VN;oi) 


flCK 


■■,Vn) 




■ ,vn;oi) 


■= EcWC7W( yi) . 


-,Vn) 




■ ,Vn;oi) 




■,Vn) 
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where the notation jj, C k denotes fij < Kj for each j = 1,...,N and c KK ^ 0. These 
results can be established by using known formulas for the action of the operators 

N B 

Ek := (2 ' 12a) 

i=l 1 

N f) 2 9 r k 8 



for k — 0, 1, 2 which for future reference we list here 



Eo^fr^ = E rT ( 2 - 13a ) 



E^ix) = \k\C^(x) (2.13b) 

1(1%-^)^ (2 - i3c) 



1 + \k\ 



2 N / 2 \ 

L> 2 Ci Q )(x) = 4C( Q) (^), 4:=-|K|(JV-l) + E;«i(«i-l--(i-l))(2.13e) 



(the action of -D can be computed from the commutator formula Do = [Eo, -Di])- Here 
the generalized binomial coefficients f *J are defined by the expansion 

CW(l+t 1 ,...,l+t J y) H ^ /«\ (?(">(*!,. 

where Cj^(l N ) has the explicit form 

= n (JV-(i-l)+a(j-l)). (2.15) 

We have also used the notation 

(note that this is the opposite of what is used in |27], [HJ but rather is that used by 

The polynomials in Q2.11| ) are referred to as generalized Hermite, Laguerre and Ja- 
cobi polynomials respectively [14]; they are uniquely specified up to normalization as the 
eigenf unctions of the operators ( |2.1| ) with an expansion in terms of Jack polynomials with 
highest weight (i.e. largest partition in reverse lexicographical ordering) c KK Cl a \ For the 
normalization we choose 

0^=2^/0^(1^, c£) = (-l)H/|«|!C7W(l JV ) and c£J = 1. (2.16) 

With this choice, for N — 1 the generalized Hermite and Laguerre polynomials exactly 
coincide with the classical Hermite and Laguerre polynomials (|2.2a|) and (2.2b) respec- 
tively, while in the N = 1 case Crfe corresponds to the Jacobi polynomial Pj. a ' b \2y — 1), 
normalized so that the coefficient of y k is unity. 
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There have been a number of studies of the generalized Jacobi polynomials G^'^ 
22], f|, |24|. In particular, it is known that these polynomials are orthogonal with respect 



to the inner product 

N 



(f\g) w --=U dy iy ?(i- yi ) b J] \yk- yj \ 2/a f(yi,...,yN)g(yi,...,yN) (2.17) 

l=l J0 l<j<k<N 

This is significant to the study of the generalized Hermite and Laguerre polynomials, as 
both are limiting cases of the Jacobi polynomials. Thus by comparing the operators 
and H {L) with and using the facts that the Jack polynomial Cfw is homogeneous of 
order |ac| and that in the expansion ( 2.14 ) the binomial coefficient is non-zero if and only 
if /i C ac [FT], we see that 



b — ci a \l N ) v 2 



'1 _ ^ In _ ^ 

>2 V 6 J '"-'2 l b 



H K {yi, ...,y N ;a) 



and 



lim 



■G^( yi /b, . . . , W&; «) = . • • , y N ] a) 



(2.18) 



(2.19) 



It thus follows that by performing the same change of variables and limiting procedure 
in (|2.17|) , we will obtain inner products for which the Hermite and Laguerre polynomials 
are orthogonal with respect to. We find that for the generalized Hermite polynomials this 
inner product is 



(f\g 



,(H) 



N 

n 

i=i 



dyie~y> n \yk-yj\ 2/a f(yi, 

l<j<k<N 



t y N )g(y x ,...,y N ) (2.20) 



while for the generalized Laguerre polynomials it is 

(f\g) {L) ■■=![ dmy?e- yi II \yk-yj\ 2/a f(y 

7 — 1 •'O \ A ^■Us' AT 



!)•••) 



VN)g(y 



i- 



yN) 



(2.21) 



l<j<k<N 

(these inner products have previously been identified by Lassalle [fZT], ^]). 



3 The generalized Hermite polynomials 

3.1 The generating function 

The starting point and key source of inspiration in our studies of the generalized Hermite 



and Laguerre polynomials is a private correspondence with M. Lassalle[|I|, in which we 
received unpublished notes containing, amongst other results, a multi-variable generaliza- 
tion of the classical generating function formula 

t^jr^ = e 2yz e- z \ (3.1) 

fc=0 

which is given by the following result. 
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Proposition 3.1 Let y :— (yi, . . . ,y^) and z :— (z\, . . . , zn). The generalized Hermite 
polynomials H K (y;a), defined in the previous section as polynomial eigenfunctions of the 
operator \2.1a\) , which have highest weight term as in ( 2.11a ) with the normalization 
specified by ( \2.1b\ ), are given by the generating function 



E^H K (y;a)Cj?\z)= rt a \2y;z)e 



P2(z) 



where 



1 C7<«>(2y)CW(s) 



N 



t l«|! ci a \l N ) 



and p 2 (z) := ^ 



(3.2a) 



(3.2b) 



3=1 



A fundamental result in Lassalle's researches is an explicit formula for the action of 
the operator (recall ( |2.13a|) ; the superscript (y) indicates operation with respect to 
the variables y) on ^ ( j a ' ) : 



N 



rt a) (2y; z) = 2 Pl (z) Q ^ a) (2y; z), where Pl (z) := £ z y 

3=1 



(3.3) 



This formula follows from ( 2.13a ) and the result [17, 20 1 

1 JL{k^ 



Pi 



(s)C«. 



X 



1 + \k\ 



E 

1=1 



K 



Cut) ( x ) 



Now in the notation of ( |2.12| ), the operator H^ H > ( |2.1aj ) is given by 

H&) = D - 2Ei 



(3.4) 



(3.5) 



Knowledge of the action of Dq on J r Q O! \2y; z) is required to prove Proposition 3.1. 
Lassalle uses the formulas (|2.13|) and ( |3.3|) to establish this action. We have observed 
that in fact the required formula can be derived from ( |3.3| ). In our derivation we make 
use of the general fact that if A^F = A&F and B^F = B&F, then 



j\iy)^(y) p = a(v)b(*) p — b(z)a(v) p — gW^W p 



(3.6) 



where the second equality follows because operators acting on different sets of variables 
always commute. 



Lemma 3.1 We have 

4 a \2y; z) = (-(N - l) Pl (z) + 2£«) .F (a) (2 2/ ; z) 
D^ 4 a \2y;z) = 4p 2 (z) ^ (2y; z) 



(3.7a) 
(3.7b) 



Proof Since = \ [E^ y) , D [ 2 y) } , using 



the Jack polynomials, and (|3.6| ) gives 



[Di*\p 1 (z)) 3r>(2y-,z) 
(-{N-l) Pl {z) + 2E^\^ a \2 y] 



the fact that is an eigenoperator for 



a 



where the second equality follows by computing the commutator. To derive the second 
result, note that = [E^\ d{% so from Q, ( pTTa]) and (§7 



[-(N-l)p 1 (z) + 2Ei z \2 Pl (z)] ^ a >(2y 1 z 



a 



A P2 (z) 4 a \2y;z). 



Let us now show how (3.7b) is used in Lassalle's derivation of ( 3.2a| ). 

Proof of Proposition 3. 1 We first want to show that H K (y; a) as defined by the generating 
function ( |3.2a| ) is an eigenfunction of ( |3.5| ) with eigenvalue — 2\k\. To do this, consider 
the action of E± on both sides of ( |3.2a| ). On the r.h.s. we have 

E{ z \rt a \2y;z)e-^ 



e - Mz)E (z) o ^) (2y . z) _ 2 p 2 (z) 4 a \2y; z)e 



■P2(z) 



= e-^E? rt a \2y; z) - ±D® ^ a) {2 y] z)e~^ 

= ~\ E VT^ J) - 2E ^) H ^ «)Cf°(*), (3-8) 

K 

where the second equality follows by using (3.7b) and noting that since E\ is an eigenop- 
erator of the Jack polynomials, the definition (3.2b) gives 

E{v\Ft\2y;z)=E[ z \Ft\2y;z), 

and the final equality follows by substituting the generating function. On the l.h.s., since 
E^ is an eigenoperator of C%*\z) with eigenvalue \k\, from the definition (3.2b) 

E A-Mv; a)CM(z) = £ <*)\k\CP(z) (3.9) 

Equating coefficients of Cjfi(z) in (|3.8| ) and ( |3.9| ) shows that H K (y; a) is an eigenfunction 
of the operator (|3.5|) with eigenvalue —2\k\ as required. 

It remains to check that H K (y; a) as given by ( |3.2a| ) has an expansion in terms of Jack 
polynomials with highest weig ht term 2^Ci a \y)/Ci a \l N ). This follows from the fact 
that to compute the coefficient of Cj^(z) in qTq (2y; z)e~ P2 ^ z \ the sum in (3.2b) can be 
restricted to partitions with modulus less than or equal to \k\. 

The strategy of the above proof leads us to generalizations of the eigenvalue equation. 
In the theory of Jack polynomials, Macdonald |25| (see also |29], |J) has given a family 
of differential operators {i}jv}j=i,...,iv which have the Jack polynomials as eigenf unctions, 
and for which the corresponding eigenvalues are known explicitly. These operators are 
given by 



D P N := f(ar l E T^u— --^ 



d d 



l A+ v 1 dz h " ' 1 dz k 



1=0 l<i 1 <i 2 <...<i 

V" ( 9 9 

x 2^ \ Zil+1 dz- "' Zip dz- 

i<i l+1 <...<i p <N U4, H+i U4 "h> 
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where A + := rii<j<fc<Af(- 2: fc ~~ z j)- Furthermore, if we define the corresponding generating 
function by 



N 

E 

k=0 



D N (X;a):=Y,X N - k D k N 



then the eigenvalues are given by 



N 



e{K,a;X) := 1[{X + N - j + 



(3.10) 



3=1 



The operator E\ is related to the Macdonald operator D\ by D l N = cxE\ + N(N — 1)/2. 
By considering the analogue of ( |3.8| ) with replaced by D-> N (j = 1, 2, . . . , N), a family 
of differential operators which have the generalized Hermite polynomials as eigenfunctions 
can be given The r.h.s. of fl3.8|) is then computed according to the Baker- Campbell- 
Hausdorff formula 



(/>v 



+ 



2! 



AT , 



Pi{z) 



,Pi{z) 



1 r 

n\ 



+ 



,?(*) 

N i 



Pi{z) 



Pi(z) / (3.11) 



Note that the sum on the r.h.s. terminates after the n-th nested commutator since the 
highest derivative in Djy has degree j. 

Following the derivation of the eigenvalue equation given in in the proof of Proposition 
3.1, and thus using (|3.3|) , ( |3.6|) and the fact that since D 3 N is an eigenoperator of the Jacks 
we have 

D^\4 a \2y;z)=D^ 4 a \2y;z) 

we can immediately deduce a family of N independent eigenoperators of the polynomials 
H K (y; a), together with the corresponding eigenvalues. 

Proposition 3.2 Let 



H 



(H) 



r>3'{y) 



+ 



1 



4 2 2! 



D, 



(id 

5 



+ 



4 n n\ 



D 



(y) 

> 



D 



(y) 

o > ' • " ) 



We have that H K (y; a) is an eigenfunction of W for each j = 1, . . . , N, with eigenvalue 



ej{n; a) given by the coefficient of X N ~ j in ( \3Jj\) . 



3.2 Consequences of the generating function 



The generating function formula (|3.2a|) can be used to deduce higher-dimensional ana- 
logues of the classical properties of the Hermite polynomials 



H n (-y) 

Ty H ^ 
2yH n (y) 



(-l) n H n (y) 
2nH n . 1 {y) 

H n+1 (y) + 2nH n ^(y). 



(3.12a) 
(3.12b) 
(3.12c) 
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Proposition 3.3 We have 



H K {-y ]a ) = {-lpH K { y] a). 



Proof Replacing y by — y in ( |3.2a| ) and using the fact that q Tq [2^] z) = Q T^\2y; fiz), 
where /i is any scalar, and that Cjf 1 ' is homogeneous of order \k\, gives 



1 



The result follows by equating coefficients of C^'(z). 
Proposition 3.4 We have 



E^H K (y;a) = 2j:( K ) 

i=l \ K (i)J 



Proof Applying to the generating function ( |3.2a| ) and using ( |3.3| ) gives 

= 2 Pl {z)Y,X-<H K { y] a)C£\z) (3.13) 



\K ! 



\K ! 



Using the formula (3.4) the final formula on the r.h.s. of ( [3.13| ) can be rewritten as 

1 ",K»< 1 " 



2E 



^(1 + N)! 



iV (i) i iv / \ 

? E( )i? K (!/;«)cS(z)=2E n E i^ (y;a)C<°>(z) (3.14) 

! t=l K « t=l 



Equating coefficients of C^^z) on the l.h.s. of fl3.13| ) and on the r.h.s. of ( |3.14j ) gives the 
stated result. 

Proposition 3.5 We have 



K 



2p 1 (y)H K (y ] a) = ajr( K{l) )^(N-i + l + aK t )H K(l) (y) + 2Y / ( " ] H (y) 



Proof From the generating function ( 3.2a ) we have 

2 E ^Pi(y)H K (y;a)CM(z) = 2 Pl {y) Q ^\2y-z)e-^ = (E^ ^ a \2y;z)) e -^ 

K \ K \- 

= (E^ Z) ^ a \2y; z)e-^) + 2 Pl (z) ^ a) (2y;z)e^ W 

= (eP + 2 Pl (z)) E |4^»(y; «)^ Q) (^) (3-15) 



n 



Using the formula ( 2.13a|) and writing 



\K ! 



as in the proof of Proposition 3.4 shows that we can rewrite the last expression on the 
r.h.s. of (|3.15|) as 



Cl a) (z) 



N 



\K\l 



i=l 



( a ) (IN 

K (i) \ l 

0) f-i N 



K 



Hi), 



H K( Jy;a)Ci a \z) 



The stated formula now follows by equating coefficients of C^\z) on the l.h.s. of (|3.15|) 
and the r.h.s. of the above equation, and using fl2.15| )) to rewrite C^(1 N ) / {1 N ) . 

Another consequence of the generating function ( |3.2a|) relates to an analogue of the 
formula (|2T 



Proposition 3.6 



We have 
1 



\n\=k 



Proof Set Z\ 



zn = c in fl3.2a| ) and note that 



Jtf\2y; c,...,c)=J2 f^ a) (2y) = e 2c ^ 



a 

K\l 



i{y) 



(the last equality follows from ( ^8|) ) to conclude 



E^ir^(y;«)ci 0) (i 



e 2c Pl (y) e -Nc 2 



E 

/c=0 



jyfc/2 k\ 



(3.16) 



The result now follows by equating coefficients of c k . 

Notice that each term on the r.h.s. of the above formula is an eigenfunction of the 
operator (3.5) with eigenvalue — 2|k|. Thus H k (^pi(y)) is also an eigenfunction of (|3.5|) 
with eigenvalue —2k. This latter fact can be checked directly, and has been observed 
previously fTI |. 



3.3 Integration formulas 



Using the generating function ( |3.2aj ) and the orthogonality of the generalized Hermite 
polynomials with respect to the inner product ( |2.20| ), a number of integration formulas 
can be obtained. In particular, we can obtain the multidimensional analogues of the 
classical formulas 

dye~y 2 (H k (y)) 2 = ^2 k k\ 



f°° dye-y 2 H k (y + x) 

J —oo 



X 



IT J-oo 

k COO 



2 r°° 2 

-= / dye~ y (x + iy) k = H k {x). 

\/7T J-oo 



(3.17a) 
(3.17b) 
(3.17c) 
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To present these analogues let us introduce the notation 



N 



dfi {H) (y) ■= II II \Vi - Vk\ 2/a dyi-.. dy N . (3.18) 

j=l l<j<k<N 

Proposition 3.7 We have 



where 



tfW := f d {H) { ) = 2 -N(N-i)/2 a7t N/2 -pr r(l + (j + l)/a) 

7(-oo,ooyv p vy; AJL r(l + l/a) 

Proof Multiplying both sides of the generating function (|3.2a|) by H K (y; a), integrating 
with respect to the measure (|3.18[ ), and using the orthogonality property of {H K (y;a)} K 
with respect to the inner product (|2.20| ) gives 

' ' C^(z) = e-^ / ,4 a) (2y;z)H K (y;a)d^ H \y). (3.19) 



\K ! J(-oo,oo) N 



Set Z\ — . . . — Zn — c, substitute fl3.16| ) in the r.h.s. of ( |3.19| ) and complete the square to 
show 



|«|! " J(-oo,oo)^ l<j<fc<V 

Now take the limit c — > oo. Since from ( |2.14|) and ( 2.11a|) 



f{H) , 

-V^ a) (c) = / n \vk - yi\ 2/a H.(y + c; «) d yi . . . d^. 

K ! ./(— oo.oo)-^ ,^.-^i.^»t 



™ d Q) (c) C& a) (l") 

the stated formula for A/"^ follows. The formula for A/q is a well known limiting case 
of Selberg's integral. 

The analogues of (3.17b) and ( |3.17c|) can be derived from the following integration 



formula. 

Proposition 3.8 We have 

j (00,00)" 



Proof Substitute the generating function ( p.2a| ) for qJ 7 ^ and integrate term-by-term 
using the orthogonality property of the {H K (y;a)} K with respect to the inner product 
( 2.20|) and the normalization integral of Proposition 3.7. The resulting series is identified 
as qJ 7 q(2z',w) according to the definition (3.2b). 



Corollary 3.1 We have 



4 a \2y; z)H K (y, a) dfi(y) = ^^ ^^ 
(-00,00)" Ck'(l N ) 
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Proof Multiply both sides of the integration formula of Proposition 3.8 by exp(— p 2 {w)) 
and substitute for exp(— p 2 (w))qJ 7 q (2y; w) using (|3.2a| ). The result follows by equating 
coefficients of Cj£*>(w) on both sides. 

Corollary 3.2 We have 

e ~" w «) = J^SW) ° F ° ' (2m - iz)C ^ >(iy) My) 

Proof By writing iw for w in Proposition 3.8, we can replace oJ r Q a \2y;w) by 
Q J-'Q a \2iy; w) , J 7 q CI \2z;w) by {2iz)w) and exp(— p%{w)) by exp(p 2 (w)). The result 
follows by using the generating function (|3.2a|) to substitute for exp(— p 2 (w)) J-'lf (2iz; w) 
and equating coefficients of Cj^'(w). 

The integration formula of Corollary 3.2 can be used to derive the analogue of the 
classical summation formula M 

~ H k (w)H k (z) ^ = J_ _ t2 yi/2 e -tZ(z*+wZ)/(l-t 2 ) e 2wzt/(l-t 2 )^ | t | < L 

fc= o k\2 k ^ v^F 



Proposition 3.9 For \t\ < 1 we have 

GW( W ,z;t) := E H K (w;a)H K (z;a) tlKl 

K Are 



1 .(i _ exp (-*— ( P2 (z) + P2 (w)) 



x CK0 



2wt z 



(1 -£2)1/2' ( 1 _ t 2)l/2 



where q = 1 + (N — l)/a. 

Proof Substituting the integral representation of Corollary 3.2 for H K (z; a) and H K (w; a) 
in the definition G^ H \w^ z;t), we see that the sum over n can be recognized in terms of 
qTq^ 1 and thus 

G (H \w,z;t) = ewW+wM 

x / d^ H \y a ) ( d^ H \y b ) 4 a) (2y a] -tw)^ {2y b] -iz) ^ a) '(2y B ; 

J(—oo,oo) N J(—oo,oo) N 

We now use Proposition 3.8 to integrate over y a . This gives 

(A/o ) 2 J (-°°,°°) N 

x / d/iW(^) ^ a) (2^; - t 2 )- 1/2 )o^ (a) (2i/ fe ; - 1 2 )- 1/2 ) 

J(-oo,ao) N 
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where the second equality follows by combining dfj,^ H \yb) and exp(t 2 p 2 (?/&)) (recall ( p.lS| )) 
and changing variables. The integration over dfi^ H \yb) can now be performed using 
Proposition 3.8, and the summation formula for G^ H \w, z\ t) results. 

Notice from (|3.16| ) that in the special case that w% — ■ ■ ■ — — c, the summation 
formula of Proposition 3.9 is entirely in terms of elementary functions: 

GW(w, z- t) = -^(1 - tY Nq/2 exp hf*l - 2tcz i + ^ ■ (3- 2 °) 

Interpretation of this result in terms of an explicit solution of the Fokker-Planck equation 
Q1.6Q with W given by (|1.3a| ) will be discussed in Section 5. 

In Corollary 3.2 a certain integral transform is applied to the Jack polynomial to obtain 
the generalized Hermite polynomial. It has been observed by Lassalle that the generalized 
Hermite polynomials can be obtained from the Jack polynomials by the action of a certain 
exponential differential operator. Thus from the formula (3.7b), we see that 



(i^ ) )V (a) (2y;.) = (p 2 (.))Vo (a) (2y; 



^4 

which after multiplication by (— l) k /k\ and summing over k gives 

exp ( - 1 -D^) 4 a \2y; z) = e^jfr^ z). 

Use of the generating function (|3.2a| ) on the r.h.s. and equating coefficients of Cj^(z) 
gives Lassalle's formula 



J^— exp ( - \D^)C^\y) = H K {y- a). (3.21) 



Comparison with the formula of Corollary 3.2, and use of the fact that {C^\y)} K forms 
a basis for symmetric analytic functions shows that for any symmetric analytic function 

/(.</)• 



7V (H) J (-00,00) 
From (|3.22|) we see that if 



/ 4 a \2y; -iz)f(iy) d^ H \y) = exp ( - -D^)f{z). (3.22) 

J (—00,00)" 4 ' 



'0 

then 



pP2(z) r 

F ^ = Tffi , , N o4 a) (2y;-iz)f(iy)d» iH Xy) (3.23a) 
A/J ' J(-°°,°°) N 

f(z) = exp(\D { z) )F(z). (3.23b) 



On the other hand, by replacing z by iz and f(ix) by F(x) we have 



/ Fl?\2y- z)F{y) d^ H \y) = exp (jD^)f(z). (3.24) 

J (— co. 00)™ N 4 ' 



M^ H) J {-00,00) 
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Comparison of (3.23b) and (|3.24 ) gives 



f(z) = — W I m o4 a \^ z)F{y) d» {H \y) (3.25) 

J\f^ n > ./ (-00,00)^ 



which is the inversion formula for the transform ( |3.23a| ) (in the case a = oo (|3.23a|) 
corresponds to the Fourier transform). 

4 The generalized Laguerre polynomials 

The generalized Laguerre polynomials, defined as the polynomial eigenfunctions of the op- 
erator (2.1b) of the form (2.11b) with normalization (|2.16| ), also satisfy higher-dimensional 
analogues of their classical counterparts. A number of these formulas have been proved in 



the case a = 2 by Muirhead |27j and for general a but iV = 2 by Yan f34 |. Below we will 
develop the theory of generalized Laguerre polynomials by presenting the analogues of the 
classical generating functions, the series expansion (2.2b), recurrence and differentiation 
formulas, integration formulas and a summation formula. In Section 6 we will identify 
the formulas known to Muirhead and Yan, as well as those which can be found in the 
work of Lassalle and Macdonald. 

4.1 Generating functions 

The classical Laguerre polynomials can be defined by either of the generating functions 



or 



{yz) a / 2 ^ o r(n + a + l) 

oo 

(i - z) ~( a+1 ) e yz ^ z -v = Y, L n(y) zn ( 4 - lb ) 

n=0 



where in ( 4.1a|) J a denotes the Bessel function. These generating functions have the 



following higher- dimensional analogues. 
Proposition 4.1 We have 



^\^\a + g; x; -z) = £ (4.2) 

[a + q\k 



where 



q := 1 + (N - l)/a, (4.3a) 

J**(a a-b b ■ x- z) - V 1 Ni^-Kl^^W^^) f4 , h) 

v T r (ox, . . . , Op, h, . . . , br, X, Z) -L^T . . . C W {1N) ^ 

with [c]W : = njLi (c - JO' ~ ■ ^ e ako ^ aue 

(n(i-^))" ( ^ ) o^° ) (-*; T 47) = E^(^«)^ ) W (4.4) 

1 z K . 
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wherej{{l-z) ■= Uf=i(l ~ Zj). 

Proof In each case we need to establish that L a K (x\a) as defined by the generating 
function is an eigenfunction of the operator (2.1b), which in terms of the notation (|2 . 1 2| ) 
reads 

H {L) =D 1 + (a+l)E -E 1 , (4.5) 

with eigenvalue and has an expansion in terms of Jack polynomials with highest 
weight term (— 1)^C^(x)/\k\\C^(1 n ). The proof of the first requirement relies on the 
identities 

(p[ x) + (a + l)E { ^) Q T[ a \a + q- x; z) = p 1 (z) F i l a \a + q- x; z) (4.6a) 
(D[ x) - E?\rt a \x-y) = (q - l) Pl (y) 4 a \x; y) (4.6b) 
with ( |4.6a| ) being established in the Appendix, and (4.6b) simply a rewrite of ( |3.7a| ) with 

y - y/2. 

First consider (14.21). We have 



E[ z) ^ a) (a + q;x; -z)e p ^ z) = e p ^ z) E{ z) f[ a) (a + q; x; -z) +p 1 (z)e p ^ z \F { 1 a) (a + q; x; -z). 

(4.7) 

Using (|4.6a ) and the fact that E^ is an eigenoperator of the Jack polynomials so that its 
action on Q Ji (x; —z) is the same as the action of E[ y \ the r.h.s. of (|4.7|) can be rewritten 

as 

_ D f) _ ( a + l) E W) ri a \a + q- x- -z)e p ^\ (4.8) 

Substituting the generating function ( [4.2]) in the l.h.s. of (|4.7| ) and computing the action of 
e[ z \ and comparing coefficients of C^(z) with (|4.8j ) after also substituting the generating 
function ( |4.2|) establishes the eigenvalue equation. The expansion of L^'(x;a) in terms 
of Jack polynomials as deduced from (|4.2j ) is given by the Proposition 4.3. Its highest 
weight term is (—iy^Cj^ > (x)/Cj^ > (l N ) as required. 

Now consider (|4.4j ). Setting yj := Zj/(1 — Zj) this is equivalent to 

Y.limc^{^-) = (U(i + y)) a+9 o^ a) (^y) (^) 

To verify the eigenvalue equation, note that E[ y ^ + E^ is an eigenoperator of C ( j° L \y/ (1 + 
y)) with eigenvalue and 

(e{ v) + 4 y) )(n(i + y)) a+q = (« + ff)pi(v)(na + y)) a+9 - ( 4 - iq ) 

Thus, applying e[ v ^ + E^ to the generating function ( |4.9| ) we have 

EM^(*;«)cJ a) (T^-) 
fi i y 

= E?+E?){j{(l + y)) a+ \4 a \-x-y) 
= (11(1 + y)) a+ \E? + E? + (a + q) P M)ort a \-x; y) 
11(1 + y)) a+q (E^ - Df* + (a + l)^))^-*; y) 
mi+v)) (4 X) - - (a + l)E^\4 a \-x ]y ) (4.11) 
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where to obtain the third equality we have used (4.6b) while the final equality follows from 
( |3.3| ). The factor involving 11(1 + 2/) i n the final equality can be commuted in front of 
the operators, since they act only on the x-variables. Use of the generating function (f4.9|) 
and comparison of the coefficient of C^\y/{1 + y)) in ( |4.11[ ) establishes the eigenvalue 
equation. The highest weight term in the Jack polynomial expansion of L®(x; a) as defined 
by fl4.4p is the same as the highest weight coefficient of Cff' (z) in the expansion of 



(na-.>r* E ^w^&> 



Mi' " era 1 , 

-(<*+?) 



Consideration of the form of the expansion of Cjf\z/ (1— z)) and (n(l — z )) in terms 
of {C^'(z)} a shows that the highest weight coefficient of Cjf'(z) in the above expression 
is (-1)\ k \C^(x)/\k\\C^\1 n ). Comparison with the coefficient of C^(z) on the r.h.s. of 
shows that L^(x; a) has the required highest weight term for its expansion in terms 



of Jack polynomials. 

In Proposition 4.1 the Laguerre polynomial is given by generating functions involving 
o^o an d o-7"i ■ It is also possible to give a generating function involving (recall 
(4.3b)) which includes both these generating functions as limiting cases. 

Proposition 4.2 We have 

11(1 - z)- c - q A a) (c + q;a + q ; -x- = J2 7~~~7^Y L a x (x; a) C[ a \z). (4.12) 



1-^ x[a + q}^ 



Proof The derivation closely follows that of ( f4.4|) above, with (4.6b) being replaced by 
the formula 

D{ x) + (a + LzlL)^ A a \c; a; x; y) = (E? + c Pl (y)) ^(c; a; x; y) (4.13) 

which is established in the Appendix. 

To derive the generating function (|4.2j ) from (|4.12|) replace z by z/c and take the limit 
c — > oo using the facts that 

J^m 1 J^{c + q;a + q; -x; - = o^t\ a + 95 ~ x \ z ), 

fim IJ(1 - z/c)- c - q = e Mz \ Urn [c + 9 ]< a) Ci a) (^c) = C< a >(z). 
The generating function (fO) follows from Q4.12|) by setting c = a. 



From the generating function (|4.2|) it is possible to deduce the higher-dimensional 
analogue of the series expansion (2.2b). 

Proposition 4.3 We have 

[a + q}i a) T ( k\ (-I)HCW(S) 
1^1' [a + ^C^) 

CfHx) = [a + g ]WC7M(F)E(-l) H f!) H! ^ ( lf (4.14b) 

ctCk 



[a + g]a 
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Proof The formula (|4.14a|) follows from the generating function formula 
ing the identity p], O 





|A|! 


CO 





by apply- 
(4.15) 



on the l.h.s. and equating coefficients of C^'{z). The formula (4.14b) follows from ([4.2|) 
by multiplying both sides by e~ Pl ^ z \ using the identity ( f4.15| ) (with z replaced by —z) on 
the r.h.s., and equating coefficients of C%*\z). 

A simple consequence of the generating function ( |4.4| ) is the analogue of the formula 
of Proposition 3.6, which is derived by following the steps of the proof of that formula. 



Proposition 4.4 We have 



■N(a+q)-l 



(pi(v)) = E L%(y;a)CM(l N ) 

\K,\=k 



Since each L®(y;a) is an eigenfunction of ( |4.5|) with eigenvalue — \k\, it follows that 
j^N{a+ q ) 1 (p 1 (y^ j s a] so an eigenfunction of (|4.5|) with eigenvalue — k (this feature can 
be checked directly). 

As our final result of this section, we note that the proof of ( |4.2| ) can be generalized to 
give a family of eigenoperators for L°(y; a), together with the corresponding eigenvalues. 
These operators are analogues of the operators of Proposition 3.2 for the generalized 
Hermite polynomials, and are derived in the same way. 

Proposition 4.5 Let 

H {L \y) ■= (d 

, (-1)" 



■)P (y) 
'n 



D? + (a + l)E^,D p N {y) 



+ 



D{ y) + (a + l)E% 



(y) 



(y) 



{a + l)E^\D p N {y) 



where D P N is the operator introduced in Section 3.1. We have that L a K (y\a) is an eigen- 
function of W L \y) for each p = 1, . . . , N, with eigenvalue e p (n; a) given by the coefficient 
ofX N ~ p in flOZU . 



4.2 Recurrence and differentiation formulas 

The classical Laguerre polynomials satisfy the recurrence relations 

xL a n (x) = (2n + a + l)Z£(a;) - (n + l)L a n+1 {x) - (n + ajL'^x) 

n 

L a r +1 ' 



m=0 

L a n (x) = L a +\x) - Lf\x) 



and the differentiation formulas 

d L a Jx) 



-L a n t\(x) 



dx 

x-^K(x) = ^-(n + a)^)- 



(4.16a) 
(4.16b) 
(4.16c) 

(4.17a) 
(4.17b) 
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The generalized Laguerre polynomials satisfy higher-dimensional analogues of these for- 
mulas. Let us first consider (4.17b). 

Proposition 4.6 We have 



E[ x) L a K (x ] a) = \K\L a K (x ] a)-^- j:( " ) U + a + —) L" 

\k\ \K (i) J V a 



Hi) 



Proof From the generating function (4.9) and the fact that e[ x ^ is an eigenoperator of 
C^\x) we have 

fj, *- > y 

= I[a + y) a+q E[ y) o4 a \-^y) 

= (E{ y) - (a + q) Pl (y/(l +y)))l[(l + y)^\^\- X - y) 
= (E? - (a + q) Pl (y/(l + y))) E LJ(x; ^\^-) 



(4.18) 



But with Zj := y 3 -/ (1 + y^), 

^c?^) = (i^-i^cf (*) 

= \^(z)- — 



E 

i=i 



(T) 




a ) 







where the second equality uses ( |2.13c| ). Substituting this expression in the r.h.s. of ( [4.18 
and using ( |3.4|) to simplify the remaining term on the r.h.s. of (|4. 18|) gives 



E?> E L;(x; oc)^^) = E L;(x; a) | MC™^) 



The result follows by equating coefficients of C^{j+ y •■ 

Next we will derive the analogue of ( [4.16a ). 
Proposition 4.7 We have 



1 JLl K 



Pl (x) L a K (x; a) = (2|«| + iV(a + q)) L a K (x; a) + — ]T ( «* + a + L« (x; a) 



i=l 



JV-l 



-(H + l)aW KW ^) (iV - * + 1 + ««;,) L" (i) (a;;a) 



Proof From the generating function (|4.9|) 



Epi(x)^(x ; a )cM(-^) = pi(x) + y)^ Q) (-^; V) 
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■4 y) + y y+\4 a) (-x; y) + A a \-^ y )E$> JI(i + y) a+q 

- + (a + q) Pl (l/(l + y))) + yf + \^ a \-x; y) 



- E$> + (a + q) Pl (l/(l + y))) £ ^(x; a)CW(^) (4.19) 

The task is now to write £?o C^(i+~) ano - Pi(V(l + ^))^t a Hl+^) as a series in 

{C/S^dy V To do this let z,- = %/(l + so that 

# = X;(l-2^ 2 ^- = 4* ) -2S? ) +4" ) and pi(l/(l+y)) = pi(l — z) = N—pi(z). 
j=i az i 

We then have 

i#)cf(_|_) = (^)-2^ + 4* ) ) C i° ) W 



f ( M c, M Q) ( ljV ) C (a)/ y )-2iuic (a) r y ) 



and 

p 1 (l/(l + y ))C(-)( T ^-) = (N- Pl (z))cf\z) 

= nc ^)(-^l_) — —-yr c ( " ) ) (^)(4.2i) 

where to obtain (|4.20| ) we have used ( |2.13a| ), (2.13b) and (|2.13c| ), while to obtain (|4.21|) we 
have used (|3.4|). Substituting ( |4.20|) and ( |4.21| ) in the r.h.s. of ( f4.19| ), equating coefficients 

of Cfio(ife) with the Lh - S - of (EUD' and use of (|2TT5| ) to rewrite Cjri(l N )/Cj$ {l N ) gives 
the stated result. 

The generalizations of (4.16b) and ( |4.16c| ) are given by the following result. 
Proposition 4.8 We have 



min(A r ,|«;|) i n 

L a -\x-a) = £ {-aY J2 p^-Al; K/a (a)L a a (x-a) (4.22) 

k/ct a vertical r-strip 
l"l I || 

L^l«{x-a) = J2a- r £ p,<f> K/ M)L a Ax-,a) (4.23) 

k/ct a horizontal r-strip 



where 



<m«)= n n S^n^wnw*))- 1 ( 4 - 24a ) 
«m«) = n i| n S^n^wnwwr 1 ( 4 - 24b ) 
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Here R K / a denotes the union of all rows which intersect n — a, C K / a denotes the union of 
all columns which intersect k — o , and h%(s), h%(s) etc. are given by (2.10b). 



Proof First consider ( 4.22| ). From the generating function (fOj) we see that 



N 

n 

3=1 



(4.25) 



Using |30 



N 



N 



-iy 



na-^) = E 1 ir-^w 



r=0 



V : 



the Pieri formula EH 



jW J(«) = r! 



E f^(«)4 a) , 

re/tr a vertical r-strip 



(4.26) 



and the relationship (|2.9|) , the r.h.s. of (|4.25| ) can be rewritten as 



N 



E E (-l) r a H k|! I^aWa)^*^ 



(4.27) 



r=0 



k/ct a vertical r-strip 



The result now follows by comparing the coefficients of Cj^(z) on the l.h.s. of (|4.25|) with 
that in fl£27|) . 

Now consider ( 4.23Q . Using the formula [ 30 1 



N 



na 



E 

r=0 



and the Pieri formula (|4.26| ) in conjugate form 



a r r\ 



n/ff a horizontal r-strip 

together with the formula (|2.9| ), we have 

M j = l <7 

OO 

= E 

r=0 « 

k/ct a horizontal r-strip 



E L:(x;a)«l ff H^^(„)C(«)(, 

At ! 



The identity Q4.23| ) follows by equating coefficients of C^ a '(z). 

The analogue of ( |4.17a| ) takes the form 
Proposition 4.9 



min(AT,|K|) 

E 

r=l 



ra 



^rip K/a (a)L a a + \x;a) 
\k ! 



k/<t a vertical r-strip 
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where ip K /o-(az) is given by (H.24dj ). 

Proof Applying Eq to the generating function (|4.4|), we have 



K j \ 1 Z / 



J 

z 



-Pi {—jW-Zi)12K + \x-,aH a \z) (4.28) 



Certainly 

Vi ZJ 3 k p^k 

If we differentiate w.r.t. t the identity 

N N (_-i\r 

n(i-^)=E^#(^ r 

i=l r=0 r - 



giving 

N 



£ 7^TTw^ ?> W* r_1 = " E ^ n(l - ht) 

r=l V fc p^k 



and set t = 1, we obtain 



Z , A (-1)' M). 



pi {yr- z ) IK 1 " *>) = E ( 4 - 29 ) 



Inserting (fl2|) back into QOp gives, after some manipulation 



N 

<x) sr^ T a (rr ,. n \ri(a) 



r=l K.CT 

k/ct a vertical r-strip 



which yields the result upon comparison of the coefficients of Cj^'(z). 
4.3 Integration formulas 

The classical Laguerre polynomial obeys integration formulas analogous to the integration 
formulas fl3.17| ) for the classical Hermite polynomial: 



k\ 



/Ve-»(L2 M ) 2 d ,= E<^±I±^ (4.30a) 

/■CO 

... . -\\^ X / y a e-\F x {a+l--xy)Lt{y)dy = x k (4.30b) 
I (a + 1) Jo 

POO 

/ y a e-%F 1 {a + l;-xy)y k dy = L k {x) (4.30c) 
Jo 



k\T(a + 1) Jo 

The higher-dimensional analogues of these formulas can be established using the gener- 
ating functions ( |4.2| ) , ( [4.4| ) in much the same way as the higher-dimensional analogues of 
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(|3 . 1 7| ) were established using the generating function (|3.2a|) . To present these results, we 
will make use of the notation 



N 



d^ {L \y) ■= n v* e_% II \vk - yj\ 2/a dyi ■ ■ ■ dy N . 

3=1 ±<j<k<N 



(4.31) 



Proposition 4.10 We have 



where 



, N-l 

A/f } := / d^ L \x) = a a-iM*-i)V«> TT 



(L.(x;a))V>(.)=^ 5 |^ 



)»/«) rr r(i + (j + i)/a)r(a + i + j/a) 



(4.32) 



T(l + l/a) 



. (4.33) 



Proof Multiplication of both sides of the generating function (|4.4p by L%(x; a) and 
integration with respect to d/i^ L '(x) gives, upon using the orthogonality of {L^}^ with 
respect to the inner product fl2.21|) , 

11(1 - z)-^ ( 4 a \-x; -r^)L a K (x; a) d^ L \x) = N^C^{z). 

J[Q,oo) N 1 — Z 

Setting Z\ — . . . — zn — c, using ( [3.161 ) and changing variables cxj/ (1 — c) =: yj this reads 



c -N(a+q) 



L e- My)/c Uy]L a J^y?) n \y k -y/ /a f[dy3=^ L) cl«\c N ). 

■ / [°.°°) JV 7=1 V C I Ki<k<N 7=1 



l<j<k<N 

The stated result follows by choosing c = 1 and noting from ( |4.14a ) that 



L a Ma) = [a + q]^/\K\\. (4.34) 
The analogues of the formulas (4.30b) and (|4.30c|) are deduced from the following 



integration formula. 
Proposition 4.11 We have 



/ o^i a) ( a + <?; x; -z a ) o^i (a) ( a + 9\ x; -z h ) d^ L \x) 

= N^ L) e-^ Za) e-^ z »\F[ a \a + q; z a , z b ). (4.35) 



Proof Substitute for eP 1 ^ Ji (a + q;x;—z s ) (s = a,b) using the generating func- 
tion ( |4.2| ) and integrate with respect to d^ L \x) term- by-term. From the orthogonality 
property of {L a K } K with respect to ( |2.21| ), only the diagonal terms in the double sum 
are non-zero, with the integral then being evaluated according to Proposition 4.10. The 
resulting sum is identified with qT^ according to the definition (4.3b). 
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Corollary 4.1 We have 



r AT {L) 
/ f[ a \a + q; x; -z a )L a K (x; a) d^ L \x) = 



-pi 



(4.36a) 



rt a \a + q; x; -4)^(1) d^ L \x) = J^ L) Ci a \l N )\K\\e-^L%(z a] a). (4.36b) 



Proof The integration formula ( |4.36a ) follows from (|4.35|) after multiplying both sides 
by e Pl ^ Zb \ using the generating function ( |4.2| ) to substitute for e Pl ^ Zb \j rt f\a + q;x; —Zb) 
and equating coefficients of Cf£*\zb) on both sides. The integration formula (4.36b) follows 
from (|4.35|) after replacing Zb by — Zb, substituting for eP l ( Zb \j~i (a+q; z a ; —Zb) using (|4.2|) 
and equating coefficients of Cjf' (zb) on both sides. 

Analogous to the Hermite case, the integration formula ( [4.35 ) can be used to derive 
the analogue of the classical summation formula valid for \t\ < 1 



5 (a + 1)» 



■Ll{z)L'{y)1» = (l-t)-~- 1 e*p 



l-t 



x+y)) F 1 [a + 1; 



xyt 



(4.37) 



Proposition 4.12 For \t\ < 1 we have 
G^(x,y;t) . ^^L jy-a) ^ 



(1 - exp ( - ^(piix) + Px(y))) rt a) (a + q; JL-j - tX 



Proof This follows by following the procedure used in the Hermite case, Proposition 
3.9. 

Notice that in the special = the above summation reduces to elementary 

functions, giving 



G iL \0,y,t) 



1 



A/n L) 



(l-t) 



-N(a+g) 



cxp 



t 



l-t 



pi(y))- 



(4.38) 



Interpretation of this result in terms of an explicit solution of the Fokker-Planck equation 
(|1.6j ) with W given by (1.3b) will be discussed in the next section. 

Finally, we note (prompted by M. Lassalle) that (|4.6a|) implies the identity 



4)M 



^(1 exp ( - flW - (a + l)E^)C^(x) = L a K (x\ a) 



(4.39) 



(the derivation parallels that of ( |3.21|) ). Analogous to ( |3.22|) , comparison with (4.36b) 
gives that for any symmetric analytic function f(x), 

77m / N A a \a + q-x--z)f(-x)d^ L \x) = D {z) - {a + l)4 z) )f(z) (4.40) 
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and hence by a similar argument as before, if 



then 



F W = 7m, ^ N o^ a \a + q;x;-z)f(-x)d^ L \x) (4.41) 



fW = -^T , ^ N oFl a \ci + q;x;z)F(x)d^ L \x) (4.42) 
(in the case a = oo ( |4.41|) corresponds to the Hankel transform). 

5 Applications 

5.1 The ground state global density 

As noted in the Introduction, the ground states of the Schrodinger operators ( |1.2| ) are, up 
to normalization, of the form e _/3H///2 where W is given by Q1.3Q - The ground state density 
in a system of N + 1 particles, Pn+i(x) say, is then given by 



< \ N + 1 



Z n +i l=1 



N 

IJ / dxi e-W*,*!,..,^) ( 51 ) 

7 1 " I 



where 

N+l 



Z N+1 ■= TT / dxi e -PW{*i,x2,.,*N+i) (5.2) 

r= i Ji 



An alternative interpretation of ( |5.1|) is as the density at the point x in the statistical 
mechanical system of N + 1 particles with potential energy W confined to the interval /, 
in equilibrium at inverse temperature (3. 

In the case W = as given by fll.3a| ), a physical argument based on the interpre- 



tation of the harmonic term as an electrostatic potential (see e.g. 0) predicts that for all 

P 

lim M p(VWx) = { Z VY ^' !f! J ! (5-3) 



N^oo V N 1 0, \x\ > 1 

This limit gives the so-called global density, and the result is known as the Wigner semi- 
circle law. 

For W = the change of variables x/j = xj gives 

PN + M = i ± ±e-^y^l[ d yi \y - Vl f e-^y^ /2 ]J \y k - y/ (5.4) 

^N+l l=1 J0 l<j<k<N 



where fj, :— a' — 1/(3 and 



N+l roc 



Zn + i:=U dme-^y^ 2 ]J \y k - y/ (5.5) 

1=1 J0 l<j<k<N+l 



The same type of electrostatics calculation used to obtain ( [5.3|) predicts that for all (3 and 
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Finally, for W = W^ J \ the change of variable sin 2 <pj = yj gives 



pn + x{v) = i^y^ii -y)^ /2 U I dm \y-yi\^' 2 {i-yi)^' 2 n 

^N+l l=1 J0 l<j<k<N 



Vh-Vjl 
(5.7) 



where \i\ := a' — 1/(3, fi 2 '■— b' — 1/(3 and 

N+l i 

z N+ i= n / dy iy ^ /2 (i- yi r^ n ii/*-v/ m) 

In this case the electrostatics calculation gives 

1 f L ^ < y < 1 

lim ^PN + i(y)=\ (5.9) 
JV-»oo TV [^0 y < 0, y > 1 

In this section we will show how for (3 even the density ( |5.1| ) in the Hermite, Laguerre 
and Jacobi cases is related to eigenstates of the operator (|1.7|) and thus the generalized 
Hermite, Laguerre and Jacobi polynomials respectively (this result is already implicit in 
earlier publications || 10, 0). Furthermore, we will show how the global density can be 
evaluated by using integral representations. 



5.2 Relationship between the density and the generalized poly- 
nomials 

Instead of considering the density directly, we proceed as in || and introduce a function 
/ depending on the auxilary variables t\, . . . , t m : 

f(h, . . . , t m ) := — I] / dy t e-f"W> JJ \\{ yi - Q ]J \y k - Vj f (5.10) 

i=i JI 8=11=1 l<j<k<N 

where the normalization is chosen so that / equals unity at the origin. For an 
appropriate choice of I and V , ( 5.10|) gives each of the densities in the Hermite, Laguerre 
and Jacobi cases for (3 even according to the formula 



Pn- 



-i(?/) = (iV+l)#^e-^/(t 1 ,...,t,)_ _ (5.11) 



Zn+i 

Let us consider each case in turn, starting with the Jacobi case. 
Jacobi case 

Kaneko KITH has shown that with 



=tp=y 



I =[0,1], =y^{l-y) x \ t := fa, . . . , t m ), \ = (3^/2, (i = l,2) 

(5.12) 

/ := /^(^l; ^2 j /3j t) as given by ( |5.10| ) is the unique solution of each of the p.d.e.'s 

f&p 2 2 BF 

t P (l - t p )— + (c' - -(m -l)-(a> + b> + l- -{m - ^ - a'b'F 



^^(Ml-M^-Ml-^-O (5.13) 

3^P 
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(p = 1, . . . , m) with 



a' = —N, m — (3, b' = ~(Ai + A 2 + m + 1) + N - 1, J = ~(Ai + m) (5.14) 

P P 

Furthermore, Kaneko (see also Yan [33]) has shown that the solution of ( |5.13 ), nor- 
malized to unity at the origin, is also given by the generalized hypergeometric function 
2 Fi 0/2 \a',b';c';t) where 

1(0/2) . . . K/3/2) 



5 Ff / 2 )( 0l , . . . , a p ; &1 , . . . , fc r ; t) := £ ^ ^ ' ' ' C<f"*\t) (5.15) 



(c.f. (4.3b) ), so that 



f^(X 1 , A 2 , 0; t) = 2 Fr z '(-N, ^(X 1 + \ 2 + m + l) + N-l;^(X 1 + m);t) (5.16) 



i09/2) 











Notice that with a' = -N, 2 F x (/3/2) (a', 6'; c'; £) as defined by (|5.15|) indeed terminates and 
gives a polynomial. To see the connection with the generalized Jacobi polynomials, we 
note that by summing the p.d.e.'s ( 5.13|) an eigenvalue equation results. The eigenoperator 
is precisely the operator (|2.1c| ) with 

N = m, a = |(Aa + 1) - 1, b = ~(A 2 + 1) - 1. 

Furthermore, from (|5.15|) and (|5.16|) has a Jack polynomial expansion of the form 
( 2.11c ) and from the definition ( 5 . 1 0|) of /, the highest weight monomial in the power 
series expansion of / is mr N my We thus have 



/«(A 1; X 2 ,(3;t) = Gp+i)/^Wi) ( ^ /2) 



(5.17) 



The tilde here denotes that the normalization in the generalized Jacobi polynomial is 
such that it equals unity at the origin. Comparison of ( |5.16|) and ( |5.17|) gives an equality 
between G and 2 iq therein. 

The formula ( |5.11| ) for pjv+i(z/) also requires the value of Qn and Zn+i- Both quantities 
are examples of the Selberg integral 



Sat(Ai, A 2 , A) := (f[ fdttttil-t^) J] 

\l=l J ° ) Kj<k 



\tk — t- 



|2A 



±<j<k<N 



JV-1 

n 

3=0 



r(Ax + i + jA)r(A 2 + i + jA)r(i + (j + i)A) 

r(Ai + A 2 + 2 + (N + j - l)A)r(l + A) 
We have 

Z N+1 = SwiPmfaPtofapft), Q N = + W/W2,/ty2). 

Substituting ( |5T7| ) and ( |5TT9|) in (|5TTI| ) gives 

Sjv+i(p>i/2, p/i 2 /2, py2) 



(5.1? 



(5.19) 



^(2/^+^-1,2/^2-1) 
XLr (Ar(3) 



(ti,...,^;/3/2) 



t 1 =...=tp=y 



(5.20) 
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Our ability to compute the global limit relies on an integral representation of 2^^^ 
(different of course from ( |5.10|) ) and thus, after equating ( |5.16| ) and ( |5.17|) , of G^ m )- This 
integral representation can be derived from the integral representation of the generalized 
hypergeometric function p3| 



2 F{ 1/X \a,\(m- 1) +Ui + l;2A(m- 1) + Vl + u 2 + 2;t) 



Sm(vi, U 2 , A) 7 [0,1] 



/ dx 1 . . .dx ml J r ^ /X \a;t;x)D u u x (x) (5-21) 
J\o,i] m 



with iJ r Q 1//A ' ) given by (4.3b), and 



3=1 



*iT n 

l<j<k<m 



Xfc Xj 



|2A 



(5.22) 



Since G in ( |5.20|) is equal to 2-^1 in ( |5.16|) with m = (3, we must set 

4 

a = -N, A = 2/P, ui = - + iii + a* 2 + N - 2, i/ 2 = -2 - ^ 2 - JV 



(5.23) 



in (ggg) . 

We note that z/ 2 is negative so that ( 5.21Q is not defined as written. However we 
can readily analytically continue the integral (|5.21 ) so that it is valid for v 2 negative 
by following the procedure detailed in [|Kj]. Thus we deform the contours [0, l] m to the 
contours C m , where C is any simple closed contour which starts at the origin and encircles 
the point x = 1 (this is first done under the assumption that v 2 is not an integer, and A 
is an integer; it is extended to all v 2 by analytic continuation and to all A by noting that 
the r.h.s. is analytic in A when it is defined, while the l.h.s. is a rational function of A in 
the case of interest (a = —N)). Furthermore, we have the formula [17 



Thus we have 

Mz/p+m-ip/p+nz-i) 



1 f 13 

— J g dx x . . .dx p Y[(l - yxi) N x 



x J-q ^ (a; t; x) 



(t a ,...,^;/3/2) 



tl,...,tjS=C 



IE 1 - cxi)' 
1=1 



(5.24) 



t 1 =...=t /3 =y 



i- xl )- 2 -^- N n w-xj 

l<j<k<P 



(5.25) 



where C is chosen so that at t = 0, G is unity. This is our desired integral representation. 
Laguerre case 

Let (|5.10|) with Laguerre weight e~^ v ^ = e _/3y / 2 y^/ 2 and integration interval I = [0, 00) 
be denoted / = p L >(fi,P; t). Comparison with the definition of f^ J > shows 

, 1 v mN 

hm f( J \(3fi/2, PL/2, P; t/L) = /^(/3/i/2, ft t) (5.26) 

Substituting (|5.16| ) and (|5.17|) for and using ( |2.19| ) and the fact that 
lim^oo 2 F[ a \a, b; c; x/b) = xFi a \a\ c; x), we thus have 



] {L \p l xl2-t) = L 



H-1+2/P, 
(Arm) 



it ] p/2) = 1 F[ m {~N^ + 2 ] t) 



(5.27) 
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where L denotes the generalized Laguerre polynomial normalized to unity at the origin. 
(The equality between and iF 1 ^ 2 ' ) has previously been given in jTIJ and the equality 
between L^i"| 2///3 and iF^^ 2 ^ has been noted in ||23|| .) Furthermore from the working in 
Toj we have 



lJ Pf /2) (-iV;/x + 2;t 1 ,...,t /3 

P 

d.T.i ■ ■ ■ H,T,a 

a J ic >) 



ti=...tp=y 



J H C ') P j=l l<j<k</3 



It, - r-l 4 /' 3 

| ^ ft ^ J | 



(5.28) 



Hermite case 

Starting with / and e~^ v ^ given by ( p.!2[ ), the Hermite case I = (—00, 00) and e - ^^ 



can be obtained by the change of variables and limiting procedure 
1- Vjs . 1/- u 



yj "^2 (1 "f ) ' *^2 (1_ L ) ' ^i/ 2 = = /5L 2 /2, L^oo (5.29) 

Hence from ( [2.18 ) and ( 5.17]) with Ai = A2 = (3L 2 /2 we see that in the Hermite case, 
is proportional to Hr N m^(ti, . . . ,t m ; (3/2) (the fact that p H ' is an eigenfunction of 
( |2.1a| ) with eigenvalue — 2N was shown in U). Thus, from (|5.1|), if we denote by H K the 
generalized Hermite polynomial normalized so that the coefficient of the highest weight 
monomial m K is unity, we have 



p N+1 (x) = (iV + l)^L e -^ 2 /2# ( tlt ... ttp .p/ 2 ) 



.ta=X 



where 



N 

n 

1=1 



d\i exp 



4e\ 2 ] n 



Afc — A.; 



P 



J I 



i=i / i<j<k<N 
N/2-Np(N-i)/4(2> K ) N /' 2 r(l + + l)/2) 



3=0 



T(l + (3/2) 



(5.30a) 



(5.30b) 



(compare (5.30b) with A/o in Proposition 3.7). To obtain a form of H^ N p^ suitable for 
asymptotic analysis, we make use of the integral representation Corollory 3.2 of H K . In 
the case of interest (k = (N 13 ), t\ = . . . — tp = x) we have 



J r o 2/P) (2y u ...,2yp; -iz h . . . , -izp) 



--z -- 



rj e - 2ixy i 
3=1 



and 



P 

■iWp) = Him] 

3=1 

so we can complete the square in the integrand of the formula of Corollary 3.2 and change 
variables to obtain 



-tp=X 



— / dm . . . dup JJ (iUj + x) N e u j Yi \ u k 

V P J]R 7=1 l<j<k<(3 



(5.31) 
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where 

V m := Af( H) (N = m,a = (3/2) (5.32) 

(It is also possible to derive ( |5.31| ) by performing the limiting procedure (|5.29|) in the 
integral representation fl5.25j ).) 

Analogous to the situation in (|5.25|) and ( |5.28| ), we note that each integration path 



along the real line can be deformed to the path C", where C" is a simple contour which 
starts at — oo and ends at oo (this is true for 2/(3 G 2> by Cauchy's theorem; it then 
follows for all values of 2/(3 that the r.h.s. is defined by noting that the r.h.s. is then 
analytic in 2/(3 while the l.h.s. is a rational function in this variable). 

5.3 The global density limit 



Using the integral representations (|5.25| ), ( |5.28| ) and (|5.31| ), the global density limits in 



(5.3), (5.6) and (|5.9|) can be computed for all (3 even. The method used in each case is 



to deform the contours so that they pass through the saddle points (for each integration 
variable there are two saddle points), and to expand the integrand in the neighbourhood 
of these points. Due to the similarities of the three calculations, we will give the details 
in the Hermite case only. 

Changing variables u\ \— > V2Nui in the integral representation ( |5.31| ), and substituting 
the result in ( |5.30a|) gives 



Zn /n Ar\('tf/V+,3tf-2l/2 .-tf/V.T! 2 



PN+1 (V2Nx) = (N + l)-^—(2NY l3N+ ^- 2 y 2 e- l3Nx 

Zn+iV/3 

X 



/ d Ul ...dupf[e- 2Nu i(i Ul + x) N J] \u k - Uj \ 4 ^ (5.33) 

P Z=l l<j<k<(3 

For each integration variable u\ the iV-dependent terms in the integrand are 



e 



-2Nv?/ iui _|_ X \N _ e -2Nuf+Nlog{iui+x) 



The exponent has a stationary point when 

= ±^(l-x 2 ) 1/2 , (5.34) 

so according to the saddle point method of asymptotic analysis we should deform each of 
the contours of integration in ( |5.33| ) to pass through u + and 

With the contours of integration so deformed, we must expand the integrand in the 
neighbourhood of the saddle points. Due to the factor Ili<j<fe</3 \ u k~ the maximum 
contribution will be obtained by expanding (3/2 integration variables (ui, . . . ,up/2 say) 
about u + and the remaining (3/2 integration variables (w/3/2+1, • • • , up) about U-. This 
specific choice is only one of the equivalent ways of dividing the integration variables 
into these two classes, so after expanding the variables with the specific choice we must 
multiply by the combinatorial factor. 

In the neighbourhood of the specified points we have 

e ~2Nuf+N \og(iui+x) 

~ exp[-2iYw! + N\og(iu± + x) - \{u - u±) 2 {AN - - - - )] (5.35) 

2 (m± + x) z 
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where on the r.h.s. u + is to be taken for j = l,...,/3/2 while «_ is to be taken for 
j = P/2 + l,...,p. Also 



Yl \uk — Uj\ 4 ^ ~ \u + — u_ 

l<j<k<P 



p n w-u 3 \^ n 

l<j<fe</3/2 /?/2+l<i<A;</3 



u fc -u/ //3 (5.36) 



Thus after substituting ( |5.35| ) and ( |5.36| ) in ( |5.33| ) we obtain 
p N+1 (V2Nx) ~ 



\J^(2N) 



Zn+iV/3 



(J3N+30-2) /2 e -/3Nx 2 -N ' f3{u 2 + +u1)+{N f3 /2) log |m++2;| 2 



U + — U- 



/3/2 



B/>/ : 



dui . . . n exp[-2iVw 2 (2 A 



N 



i=i 



2{iu + + x) 2 



)] n l«* 

l<i<fc</3/2 



it 



I4//3 



(5.37) 



To simplify ( p.37| ) note that a simple change of variables gives that the last line is 
equal to 



1 



p-i 



\2N 



N 



(Vp/2f 



2(iu + +x) 2 I 

where Vp/2 is defined by (|5.32j ). Now suppose x < 1 so that u*_ = —u + . Using (|5.34 ) we 
then have 



2N- 



N 



2 (mi i + x) 



4iV(l-x 2 ) 1 / 2 , 



Ux — U- 



(1-* 2 ) 1/2 , 



1 1 

2,2 2 I • , I 1 

U + + U_ — - — X , \IU + + X\ — -. 



Making these substitutions in (|5.37|) shows that 

p N+1 (y/2Nx) ~ (1 - x 2 ) 1 '^ ?" (V ff )2 f e -W/*2-WWW N fW+W (5.38) 

Zn+i vp 

To simplify the x-independent terms in ( |5 . 38 ) we note from the specific formula (5.30b) 
and Stirling's formula that 



(3 



(i+PN)/2^ 27l yi/z 



r{i + p/2) 



T(l + (N + l)/3/2) 



IW2 + 1) 



7T 



2 PN/2-l/2 N -(p/2)(N+l)-l/2/p i 2 y{P/2) e Np/2 (5.39) 



Also, from (|5.32|) and straightforward manipulation of the explicit formula in Proposition 
3.7 we have 



Vp ~ w } r(i + /3) 

Substituting (|5.39| ) and (|5.40|) in (|5.38|) gives 

\/2N 

p N+1 (V2Nx) ~ (1 - x 2 ) 1/2 , \x\ < 1 

7T 



(5.40) 



(5.41) 
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which is precisely the formula ( |5.3|) for |a;| < 1. 

For the intervals \x\ > 1, instead of repeating the working of the expansion about the 
saddle points (which are both pure imaginary in this case), we note from the result (|5.41 
that 

2N r 1 



p(V2Nx) d(V2Nx) 



71 



'l-x 2 ) 1/2 dx = N 



But from the definition of the density it is non-negative and satisfies the normalization 

p(V2Nx)d(V2Nx) = N. 



Thus we must have 



as predicted by (K 



p{V2Nx) 
V2N 



0, for Ixl > 1, 



5.4 Initial value problems 

The summations G^ H \w, z;t) in Proposition 3.9 and G^(x,y;t) in Proposition 4.12 are 
essentially the Green functions for the solution of the Fokker-Planck equation ( |1.6| ) with W 
given by (|1.3a|) and (1.3b) respectively. To see this, we first recall that P = G(x^ \x; r) is 
the Green function solution of (P..6I) if it is the solution which satisfies the initial condition 

P(x: t, 



N 

i=i 



By appluing the transformation (|1.4| ) the Fokker-Planck equation can be written as 
the Schrodinger equation (|1.5|) , where t = r/i/3. In general the Green function solution of 
the Schrodinger equation, G(x < -°- ) |x; t) say, may be written in terms of the eigenvalues and 
eigenf unctions of H. Thus suppose {ip K } K is a complete set of orthogonal eigenfunctions 
of H with corresponding eigenvalues {E K } K . Then the method of separation of variables 
gives 

G(x^t)=^ Mxi y K{X) e-« E * 
Thus for the Fokker-Planck equation 

G{x^\x- r) = e TEo ^-^^G(x^\x; r/i(3) (5.42) 

Now, for the Schrodinger operators ( |1.2a| ), (1.2b) 

ij)W(x) = ip { H) (x)HJx/^;a) E[ H } = E ( H) + - \k\ (5.43a) 

a 

= #^)4 a ' /a ~ 1/2 V/a; «) eP = E& L) + -\k\ (5.43b) 



a 



where 



N 

4 H \x) := n^ 2/2Q n i**-*;i i/o 

3=1 l<j<k<N 
N 

4 L \x) ■= n< /a ^ 2/2a n \4-^\ i/a 

3=1 ^<3<k<N 



3 I 
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Substituting ( |5.43a| ) and (5.43b) into (|5.42|) and comparing with the definitions of 
z\ t) and y; t) shows that 



& h \xW\x;t) 



a -N q /2 (^), 



,r)| 



(0) 



a,x/\/a;e T ) 



(5.44a) 



x 1 - -* la;; r) 



^(aVa-l/2+9) ^ ^ G (V (( x (0) ) 2 /a> % 2 / Q . e ~2 T) \ (5.44b) 



From (|3.20|) and ( |4.38| ) we see from (|5.44|) that for some initial conditions it is possible 
to express and in terms of elementary functions. Thus for x^ = c (i.e. xf* = 



x 



(0) 
N 



c ) in the Hermite case, from (|3.20|) and ( 5.44a|) we have 



A0)-_ 



-2t- 



A/f> 



-7Vg/2 



x exp 



1 



N 



all — e 



n 



|2/a 



(5.45a) 



j=l / l<j<k<N 

while for x*- -* = in the Laguerre case, ( |4.38| ) and (5.44b) give 



& L \x^\x;r) 



N 

X []l 

3=1 



(0) =0 N {L) 

exp 



[ail — e 

a=a'/a-l/2 V v 



2r^ 



2a' /a 
j 



'/a-1/ 
1 



-N(a'/a-l/2+q) 



N 



2r\ 51 



X , 



n 



V C ) 3=1 ) l<j<k<N 



\X k X j 



2|2/q 



(5.45b) 



Now that these explicit solutions have been revealed, they can be verified independent of 
the theory of generalized classical polynomials, by direct substitution into (|1.6|) with the 
appropriate W . 

Another consequence of (|5.44j) is that it implies the asymptotic small-r behaviour of 



( )/X (0) Xx ( \ (x (0) ) 2 x 2 

O^^S^TiJ and ^(a + qr;-—; — 

Thus in general, as r — > the asymptotic solution of the Schrodinger equation (|1.5|) (with 
t = r/i/3) is given by 



G(x (0) \x; r/i/3) 



P \N/2 



N 



AlIT 



lie 

3=1 



Substituting this in ( |5.42|) , substituting the result in (|5.44j ) and using Propositions 3.9 
and 4.12 gives 
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-N/22N(N-l)/2aj\J-( H ) 



| :r (0; _ x \ 

^Hl<j<k<N( x j ~ x k){ x j — x k )/t) j 



N 



(0) / 



— a n e ■ 
i 



(5.46) 



and 



F[ a) [a + q; 



(x^) 2 x 2 ' 



7T 



-N/22N{a+l/2)+N(N-l)/aj\f(L) 



2t ' 2r 



(ni<, 



<k<N\ X j 
N 



x 2 )((x 



(0)^2 



(o)> 



>)/r 



l/a 



TV 



x n(w7T)" (a+i/a) n e<B ^ 



(0) / 



(5.47) 
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where it is assumed X\ < ■ ■ ■ < xn and xf 1 ^ < ■ ■ ■ < x$ . In the case a = 2 these 
asymptotic formulas are known in the mathematical statistics literature (see e.g. p6| ) 

6 A brief literature survey 

To our knowledge, the generalized classical polynomials were first introduced by Herz 



15 1 in the case a = 2 via integral formulas over measures associated with spaces of or- 
thogonal matrices (however it should be noted that what Herz calls generalized Hermite 
polynomials do not correspond to the generalized Hermite polynomial we have consid- 
ered). Constantine and Muirhead extended the work of Herz on the generalized Laguerre 
polynomials in the case a = 2, and derived the formulas (|4.2|) [f27j , ex. 7.19], (|4.4j) |27| , 



Thm. 7.6.3], QQ2|) |27|, ex. 7.20], (plaQ |27|, eq. 7.6(4)], QO^ ) 0, Thm. 7.6.5], QOBD 



27 , Thm. 7.6.4] and Proposition 4.12 [[27], ex. 7.21] in that case (in comparing formulas 
it should be noted that Muirhead adopts the normalization c KK = (—1)^/C^(1 N ) which 
is \k\\ times the normalization we have used). For general a and N = 2 Yan |34| derived 
(0) @, eq. (5.6)] flOIaD @, last eqn. p. 251], ( ggg ) @, eq. (5.11)] and (4.36b) @, 
eq. (5.13)] (Yan uses the same normalization as Muirhead). For general a Lassalle |19| 
has reported the results ( |4.14| ) and ( |4.32| ), and simultaneous to our investigations has ob- 
tained the results ( |4.2| ) and ( |4.39|) - Q4.42|) (Lassalle uses the normalization L"(0; a) = 1). 



For values of a corresponding to Jordan algebras, generalized Laguerre polynomials have 
been studied by Dib ||. 



In an unpublished handwritten manuscript Macdonald |24] has derived some prop- 
erties of the generalized classical polynomials. His results for the generalized Laguerre 
polynomials, which overlap with the same equations of ours as does the work of Yan, are 
typically proved for a = 1/2, 1 and 2, and are conjectured to remain valid for general a. 
The validity of a number of the results in p3 for general a rely on a conjecture for the 



so called generalized Laplace transform of the Jack polynomial 

N 

'[0,oo) 



/ ^ a \-x ]y )C^\x)~[[x a 3 J] \x k -x j \V° l dx 1 ...dx N 

• /f0 ' oo)JV j=l l<j<k<N 

N i 

= [a + q t ] c^(i N )Y[y7 {a+9) ci a) (-) (6.1) 



This conjecture can be proved using results contained herein. First we calculate the 
generalized Laplace transform of the generalized Laguerre polynomial: 



N 

TT I . . . 1 2 ■''<! 7 . 1 " \ 



r N 

/ o^ r o X \-x;y)Ll(x;a) TT TT \x k - x j \ 2/a dx 1 . . .dxj 



i=i i<j<k<N 

N 



= ^ L) I[y7 {a+q) cl a \i--), (6.2) 

3=1 y 

which follows from the first equation of the proof of Proposition 4.10 after noting from 
flyj) that 



z , 1 



A — 1 -L 



rr.. \ — 

o 



J~n ( X', 



* z 1-2 
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and writing 1/(1 — 2;) =: y. We now use (4.14b) and multiply ( |6.1|) by a suitable a- 
dependent factor so that after summing over a we can replace L a a {x\a) on the l.h.s. by 
Cf'(4 On the r.h.s. we then have 



[a + <^ Q 



Substituting the value of from (|4.32| ) and using Q2.14j ) to compute the sum gives 
( |6.1| ) as required. 



The generalized Hermite polynomials of the type considered in this paper appear to 
have been first considered by James |L|] in the case a = 2. Subsequently, for general a 



Lassalle |T| noted the orthogonality with respect to the measure (|2.20|) , the normalization 
of Proposition 3.7 and the property of Proposition 3.3. Furthermore, in handwritten notes 
Lassalle |19|] has established Proposition 3.1 and has stated Corollaries 3.1 and 3.2, ( p. 21 



( |3.22j ) and ( |3.25| ). Also given in the notes is an explicit formula for the coefficients cffi 



in ( p.lla| ). Macdonald [24j] has also considered properties of the generalized Hermite 



polynomials in the form of conjectures based on derivations in the cases a = 1/2, 1 and 
2. He has obtained the normalization of Proposition 3.7, the property of Proposition 3.3 
and the integration formula of Proposition 3.8 and the generating function of Proposition 
3.1. 

M. Lassalle has pointed out to us that the exponential operator formulas (|3.21|) and 



( |4.39| ) imply an intimate connection between the theory of generalized Hermite and La- 



guerre polynomials and theory developed by Dunkl [y, [7] • Inspection of these works show 
that this is indeed so. The Hermite case is the most straightforward, which in the language 
of @, [7|] corresponds to the root system A N . Dunkl introduces the operators 

d 1*1- Ma 

where Mij is the operator which exchanges coordinates X{ and Xj. When acting on func- 
tions symmetric in x%, . . . ,xn these operators are related to D (recall ( p,12| )) by 

N 
i=l 

Also introduced is the pairing [p, q] H . For polynomials p and q homogeneous of the same 
degree 

\p,q] H := P {T x )q{x), (6.4) 

where p(T x ) means that each variable Xi in p is replaced by Tj (the ordering within the 
monomials does not matter since the operators {Tj} commute), while if the degrees differ 

[P, q]h ■= 0. 



This pairing is intimately related to the exponential operator in ( |3.21| ). Thus, as noted in 



[19], it follows from 0, Thm. 3.10] that for homogeneous symmetric polynomials p and q 



x . 
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From ( p.21|) and the orthogonality of {H K } with respect to the inner product ( 2.20|) we 
immediately have the result that the Jack polynomials are orthogonal with respect to the 
pairing fl6.4j) : 

where we have used Proposition 3.7 and ( j2.9| ) 

Dunkl also introduces a kernel K(x,y), which for the root system An and p a sym- 
metric homogeneous polynomial has the property @, Prop. 2.1] 

(we have changed variables x h- > y2x, y i— > \[2y and replaced K(\/2x, \/2y) by K(x,y)). 
Comparison with Corollary 3.1 (after substituting for H K using (|3.21|) ) gives the explicit 
formula 

K(x,y)=^ a) (2x-y). 
In the Laguerre case there are analogous connections with the work of Dunkl, with the 



underlying root system now being B^. The operators Tj are now given by (see e.g. |T3 



T, = #- + i E (i-gg + 1 - S ' S ' M ") + i±iZg(l - 5.) (6.5) 

I../ . ^ ■ ' _■ ^ JU i J i I j 

where the action of the operator Si is to replace the variable Xi by —Xi. When acting on 
a function / symmetric and even in xi, . . . , Xn these operators are such that 

N 



1=1 



Using (|6.5| ) and setting P(x) := p(x 2 ) and Q(x) = q(x 2 ), || Thm. 3.10] gives 

[ P> Q] L = _L f ( e -(D 1+ (a+l)E )\ ( e -(^Ha + l)E ) q \ dfl (L)( x ^ 



where is defined by ( |6.4| ) with Tj specified by (|6.5|). Orthogonality of {L K } with 
respect to (2.21) and use of ( |4.39|) then gives 

[Jl a \ 4 a) U = «-' K| J.[a + #4^)^- 

Also, (|4.36a|) with the substitution of (|4.39| ) and the change of variables 
gives the kernel K(x, , y) in the B N case as JF^^(a + q; x 2 ; —y 2 ). 

In the context of the Calogero-Sutherland model the expansion of the generalized 
Hermite polynomials in terms of monomial symmetric functions has been considered by 
Ujino and Wadati |H|], and a Rodrigues-type formula has been obtained [[3^], analogous to 
that recently given by Lapointe and Vinet [|18[ for Jack polynomials. Also Polychronakos 
2~8f has recently considered the monomial expansion of Hnk\(x; a) and given its normal- 



ization. To our knowledge there have been no previous discussions of the generalized 
Laguerre polynomials in the context of the Calogero-Sutherland model. 
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Regarding the global limit of the density computed in Section 5, we know of no other 
works which consider this limit directly. However, using techniques from potential theory 



Johansson [15] has recently proved that, in the Hermite case, for all (3 > 




2Nx)p(V2Nx) dx = - / f(x)Vl - x 2 dx. 
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(for (3 = 2 this result was first given in || using a mean-field approach) for any continuous, 
bounded /. The analogous result in the Jacobi case has also been obtained by Johansson 



| I6| . These results establish that the smoothed density is that predicted by electrostatics, 



while our result establishes pointwise convergence to the electrostatic prediction. 
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Appendix 

In this appendix, the equations ( |4.6a| ), (4.6b) and ( (4.13[ ) are derived as special cases of a 
p.d.e. satisfied by 2^ r { a \o',b;c;x;y) (recall (4.3b)). 

Proposition A.l 

Let 2 T[ a) (a, 6; c; x;y) be defined by (4-3b). This function satisfies the p.d.e. 

D[ x) F+(c- 4 X) F-(a + b- F - ^ F = ab Pl (y) F, (A.l) 



where D^, Ej. are defined by and r] 2 '■— \[F>2, E 2 \, and is in fact the unique solution 

of the equation of the form 

C W(x)CW (y) 



Proof We follow the method of Constantine and Muirhead|2^] in the case a = 2. With 
F given by (|A.2j) , from (|2.13|) and ( |3~4"D we have 



D ,, F = Eg(^)(«, + ^)^C<»^ (A.3a) 

N 



EE^D^Sy^fe)^ (A-3b) 



a J CM(1 N ) 

Pl(y)F = TTR?S( K r) m$) c "' <»> a- (A-3.) 
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Substituting ( |A.3|) in ( |A.1|) , equating coefficients of C^\x)/C^\l N ) and then equating 
coefficients of C$ } (y) ( K ") gives 



C + ^ - - ) A K (i) 

a 



\K\ 



a + Ki 



i - 1 



a 



b + ^ 



i - 1 



a 



A K 



This is a first order difference equation and so has a unique solution once the initial 
condition (Aq = 1) is specified. It is straightforward to verify that the solution is 



A K 



[4 C 



(«) 



The equation (|4.13|) for iT^ \a\ c; x; y) follows from ( |A.1[) by changing variables y — > 
and then taking 6 — > oo. The equation (4.6b) for J r Q a \x;y) follows from that for 
^^(a; c; x; y) by setting a — c — (N — l)/a, while the equation ( |4.6a| ) follows from 



( |4.13| ) (with a and c interchanged) for iJF^^a; c; x; y) by changing variables y — >• y/a and 
taking a — >• oo. 
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